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Chapter 4 ® 
The Fundamental of Electrodynamics \parackts 


Electromagnetism is one of the four fundamental interactions of nature, along with 
strong interactions, weak interactions, and gravitation. It is the force that causes 
electrically charged particles to interact; the areas in which this happens are called 
electromagnetic fields. Electromagnetism is the force responsible for practically all 
the phenomena encountered in daily life (with the exception of gravity). Ordinary 
matter takes its form as a result of intermolecular forces in matter. 

Electromagnetism is also the force that holds electrons and protons together 
inside atoms, which are the building blocks of molecules. This governs the processes 
involved in chemistry, which arise from interactions between the electrons orbiting 
atoms. Electromagnetic field theory has been and will continue to be one of the most 
important fundamental courses of the electrical engineering curriculum. It is one of 
the best-established general theories that provides explanations and solutions to 
intricate electrical engineering problems when other theories are no longer 
applicable. 


4.1 Introduction 


From our knowledge of college electromagnetics, we have learned that the integral 
form of Ampére’s Law are the magnetic fields because of a current distribution that 
satisfies the following relationship (Fig. 4.1): 


| A-di= | J -fida (4.1) 
Cc s 


The second kind of field that enters into the study of electricity and magnetism are 
magnetic fields or, more appropriately, the effects of such a field. It has been known 
to humans since ancient times, when the effects on the naturally occurring, perma- 
nent magnetic (Fe30,4) were first observed. Magnetic fields even were used for 
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Fig. 4.1 Electricity- 
magnetl ~Ceoes009>. 
CC 888888855 


medical treatment purposes in the early history of Egyptian civilization. The dis- 
covery of the orientation of North and South while seeking properties of this material 
had a tremendous influence on early navigation and exploration of the world via both 
sea and land. 

Early studies done by Oersted led to his discovery of the phenomenon that an 
electric current does produce a magnetic field; this allows us to have a better 
understanding of applications of magnetism. That work, together with later efforts 
and research by Gauss, Henry, Faraday, and others, brought the magnetic field into 
better focus and a partnership with the electric field. The subject of this book is based 
on such a foundation. 











4.2 Maxwell’s Equations and the Electric Field 
of the Electromagnetic Wave 


Although Maxwell’s equations were formulated over 100 years ago, the subject of 
electromagnetism has not remained static. Much of modern physics (and engineer- 
ing) involves time-dependent electromagnetic (EM) fields in which Maxwell’s 
displacement current plays a crucial role. Maxwell’s equations contain all the 
information necessary to characterize EM fields at any point in a medium. 

To understand the behavior of materials’ response to a high-power laser or radar 
beam one needs to consider the coupling of the laser energy or, in our case, the EM 
wave with materials. Therefore, we first need to know the optical reflectivity, R, and 
the transmissivity, 7, for light incident on a surface that divides two semi-infinite 
media. To understand reflectivity, we must use some general results from the theory 
of EM waves. 

For EM fields to exist they must satisfy the four Maxwell equations at the source 
where they are generated, at any point in a medium through which they propagate, 
and at the load where they are received or absorbed. Because the fields must satisfy 
the four coupled Maxwell equations involving four unknown variables, we first 
obtain an equation in terms of one unknown variable. Similar equations can then be 
obtained for other variables. We refer to these equations as general wave equations. 
It can be shown that the fields generated by time-varying sources propagate as 
spherical waves. In a small region far from the radiating source, however, the 
spherical wave may be approximated as a plane wave—that is, one in which all 
the field quantities are in a plane normal to the direction of its propagation (i.e., the 
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transverse plane). Consequently, a plane wave does not have any field component in 
its direction of propagation (i.e., the longitudinal direction). 

First, we seek the solution of a plane wave in an unbounded dielectric medium 
and show that the wave travels with the speed of light in free space. We then consider 
the general case of a finitely conducting medium. We show that the wave attenuates 
because of a loss in energy as it travels in the conducting medium. Finally, we 
introduce the concept of reflection and transmission of a plane wave when it leaves 
one medium and enters another. 


4.3 Wave Equations for Electric and Magnetic Fields 


In the regions of space where no charge or current exists, Maxwell’s equations read: 


V-E=0 (i) 

V-B=0 (ii) 

Vx E= oe (iii) (4.2) 
V x B= p€0 (3) (iv) 


where 


E = the electric field 

B = the magnetic field 

Mo = the constant Biot-Savart’s Law, known as permittivity of free space: 42 x 10’ 
N/A? 

€ 9 = the constant Coulomb’s Law, known as permittivity of free space: 8.85 x 107 
C?/(N - m’) 


12 


The preceding equations constitute a set of coupled, first-order, partial differential 
equations for EF and B. They can be decoupled by applying curl to (iii) and (iv): 


f OE oe 
Shay es) Ho€o On 
ox (Fx B) = 9(9-B)_V OF 
x (V x B) = V(V-B)-—V x | woe By 
: (4.4) 
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=> 


Or, because V- E=Oand V- B=—O, 


2 = 
VWE= —— VWB= os 4.5 
= Ho€&0 Or = Hoo Or ( : ) 


Equation 4.5 is a demonstration of separations between E and B, but they are 
second order [1]. In a vacuum, then, each Cartesian component of EF and B satisfies 
the three-dimensional wave, as follows: 


A oF 


2 $a eee 
a ee or 


(4.6) 

This supports standard wave equations in general within Cartesian form, which is 
like a classical wave equation of small disturbance on the string, where v represents 
the speed of propagation and is given by [1]: 


v=4/— (4.7) 


where y is the mass per unit length. This equation permits as solutions all functions 
of the form in Eq. 4.8, for example, in z-direction of propagation [1]: 


f(z t) = g(z—0t) (4.8) 


This mathematical derivation can be somewhat depicted as in Fig. 4.2 at two 
different two times, once at f = 0 and later at time t. Each point on the wave form 
simply shifts to the right by an amount vf, where v is the velocity. Then we can see 
that Maxwell’s equations imply that empty space supports the propagation of EM 
waves traveling at a speed 


1 


= 3.00 x 10° m/s = 3.00 x 10°m/s =c (4.9) 
V Foo 





which is precisely the velocity of light c [1]. Maxwell himself was astonished by this, 
and he noted that a “wave can scarcely avoid the interference that light consists in the 
transverse undulations of the same medium which is the cause of electric and 
magnetic phenomena.” 


Fig. 4.2, Wave propagation 
in z-direction [1] 





Ny 
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4.4 Sinusoidal Waves 


Of all possible wave forms, the sinusoidal one as defined in Eq. 4.10 is, for good 
reason, the most familiar. Figure 4.3 shows this function at time t = 0 with A as the 
amplitude of the wave. 


f(z, t) =Acos [k(z — vt) + 4 (4.10) 


(This is positive and represents the maximum displacement from equilibrium.) The 
argument of the cosine is called the phase, and 6 is the phase constant, one can add 
any integer multiple 27 to 6 without changing f(z, t). Ordinarily, we use a value in the 
range 0 < 6 <2z. A point known as the central maximum on the curve takes place at 
z = ot — O/k when the phase is zero. If 6 = 0, the central maximum passes the origin 
at time tf = 0; more generally, 6/k is the distance by which the central maximum (and 
therefore the entire wave) is “delayed.” Finally, k is the wave number. It is related to 
the wavelength A by Eq. 4.11, for when z advances by 2z/k and the cosine executes 
one complete cycle: 


_ Qn 


A 
k 


(4.11) 
As time passes, the entire wave train proceeds to the right at speed v. At any fixed 
point z, the string vibrates up and down, undergoing one full cycle in a period, that is, 


rat 


=; (4.12) 


If we now introduce frequency v as the number of oscillations per unit time and show 
it in the form of Eq. 4.13, then we have: 

1 kv ov 

v=s=2—=- 4.13 
T 2am vy ( ) 
For our purposes it is better to write Eq. 4.10 in a refined form and present that in 

terms of angular frequency, w, given in the analogous case of uniform circular 
motion that represents the number of radians swept out per unit time: 


@ = 2av= kv (4.14) 
Fig. 4.3 A sinusoidal wave Central 
propagation [1] maximum 
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Fig. 4.4 Sinusoidal wave {z,0) Central 
traveling to the left [1] PA maximum 


ie : 


Then the new form of the sinusoidal wave in terms of @, rather than v, is given as: 








f(z, t) = Acos (kz — wt + 6) (4.15) 
If both k and @ travel to the left, then Eq. 4.15 can be written in this form [1]: 
f(z, t) =Acos (kz + wt — 6) (4.16) 


This is in correspondence to and consistent with the previous convention that 6/k 
will represent the distance by which a wave is “delayed” because the wave is now 
moving to the /eft, and delay means a shift to the right. At t = 0, the wave looks like 
that in Fig. 4.4, and because the cosine is an even function, we can just as well write 
Eq. 4.16 in the form of Eq. 4.17: 


f(z, t) = Acos (—kz — wt + 6) (4.17) 


Comparing Eqs. 4.16 and 4.17, we conclude that we could simply switch the sign of k 
to produce a wave with the same amplitude, phase constant, frequency, and wave- 
length traveling in the opposite direction. 

From what we have learned in our complex variable course in college, about 
complex in view of Euler’s formula, we have: 


e — cos@+ ising (4.18) 
We can now write Eq. 4.17 as follows: 
f(z,t)Re [Actor (4.19) 


where Re[&] denotes the real part of any complex number, such as €. This allows for 
the introduction of the complex wave function in Eq. 4.20: 


f (z,t) = Ae) (4.20) 


With insertion of the complex amplitude, A = Ae® absorbs the phase constant. The 
actual wave function is the real part of /: 


f(z,t) = Re[f(z,1)] (4.21) 


Knowing f , it is a simple matter to find f- 
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Example 4.1 Combine two sinusoidal waves of f; and fo. 


Solution Let us write the following function fy as 


fs =f, +f2 = Re(f,) +Re(f,) = Re(f + f2) = Re(fs). 
with f 3= f it rg 2. You may simply add the corresponding complex wave 


functions, and then take the real part. In particular, if they have the same frequency 
and wave number, 


fs =n Aye”) zig Ane= 08 = Agee?) 
where 
A3 = A\ iy or A3e" = Ae + Are 
Now we are going to figure out what A; and 63 are as follows: 
(A3)” = (A3e) (A3e~) = (Aye! + Ape!®)(Ae! 4+ Are®) 
se (Ai)? Je (A2)” + AjA2 (ee ai ei e'2) 
= (A) 
A3 = (Av? + (Ay)* + A,A> COs (5; aa 62) 
Az = Are + Aye 











(Az)? + A}A2 cos (5; — 62) 





A3(cos 63 + isind3) = A2(cos 6 + isin d2) + Ai(cos 6; + isin 6)) 





= (A; cos 6; + Az cos 62) + i(A; sin 6; + Az sin 62) 
__ A3sin (63) A; sin (6,) + Ap sin (62) 
~ A3cos (53) Ay cos (6;) + Az cos (52) 


5 asi Aj sind, + A> sin 62 
= tan 
; A, cos 6; + Ao cos 6) 


tan (63) 














As we can see the combined wave still has the same frequency and wavelength is 
given by 


f3(z, t) = A3 cos (kz — wt + 63) 


4.5 Polarization of the Wave 


Polarization, also called wave polarization, is an expression of the orientation of the 
lines of electric flux in an EM field. Polarization can be constant—that is, existing in 
a particular orientation at all times, or it can rotate with each wave cycle. 
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Vertically 
polarized 
output 


Vertical filter 


Horizontally 
polarized 
output 


Horizontal filter 


Little or no output 


Vertical filter and horizontal filter 


Fig. 4.5 Polarization wave 


Polarization affects the propagation of EM fields at infrared (IR), visible, ultra- 
violet (UV), and even X-ray wavelengths. In ordinary visible light there are numer- 
ous wave components at random polarization angles. When such light is passed 
through a special filter, the filter blocks all light except those having a certain 
polarization. When two polarizing filters are placed so a ray of light passes through 
them both, the amount of light transmitted depends on the angle of the polarizing 
filters with respect to each other. The most light is transmitted when the two filters 
are positioned so that they polarize light in the same direction. The least light is 
transmitted when the filters are oriented at right angles to each other (Fig. 4.5) and in 
vector presentation, as in Fig. 4.6. 

Electromagnetic waves are considered transverse waves in which the vibrations 
of electric and magnetic fields are perpendicular to each other and to the direction of 
propagation (Fig. 4.7). These two fields change with time and space in a sinusoid 
fashion. Generally, only the electric field is represented, related to the propagation 
direction, because it is with the electric field in which detectors interact—for 
example, eyes, photographic film, charged coupled device (CCD) . 

Visible light makes up just a small part of the full EM spectrum. Electromagnetic 
waves with shorter wavelengths and the higher energies include ultraviolet light, 
X-rays, and gamma rays. Electromagnetic waves with longer wavelengths and lower 
energies include infrared light, microwaves, and radio and television waves. 
Table 4.1 is a presentation of types of radiation along with their wavelength range. 

The polarization of an EM wave refers to the orientation of its electric field, E. 
When the direction of Eis randomly varying with time on a very fast scale, related to 
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Xx, xX 
Vv Vv 
z. Zz 
yw ; 


Horizontal polarization 


Vertical polarization 





Polarization vector 


Fig. 4.6 Vector form of polarization wave [1] 
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Fig. 4.7 Electromagnetic wave propagation 


Table 4.1 Type of radiation 
and wavelength range 

















Type of radiation Wavelength range 
Gamma-rays <l pm 

X-rays 1 nm-1 pm 
Ultraviolet 400 nm-1 nm 
Visible 750 nm—400 nm 
Infrared 2.5 mm—750 nm 
Microwaves 1 mm-25 mm 





Radio waves 





>l mm 





the direction of propagation, the wave is considered non-polarized. In case of a 
linearly polarized wave, the electric vector has a fixed orientation related to the 
propagation direction as shown in Fig. 4.8. The polarization of EM waves can be 
produced by absorption, scattering, reflection, and birefringence. 
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Fig. 4.8 Electric field 
orientation for polarized and non polarised wave 
non-polarized EM waves 


polarised wave 


A linear polarizer is a device that only allows electric field components parallel to 
a certain direction (called the polarization axis) to pass through. Any EM wave that 
comes through such a polarizer is polarized in the direction of the polarization axis. 
After leaving the polarizer, the wave’s polarization (i.e., E-field direction) does not 
change, and the wave is thought to be linearly polarized. If the linearly polarized 
light passes through a second polarizer, the transmitted intensity, 1(@), of the waves 
as it leaves the second polarizer is given by the Malus Law. 

The specific manner in which a beam of an EM wave at a specific polarization is 
reflected (and refracted), at an interface between two different media, can be used to 
determine the refractive index of the solid. Specifically, for a particular interface 
there is a specific angle of incidence (relative to the normal vector of the surface), 
called the Brewster angle, which is related to the refractive index of a material. At 
this angle the reflection coefficient of light polarized parallel to the plane of inci- 
dence is zero. Thus, if the incident light is non-polarized and impinges on the 
material at the Brewster angle, the light reflected from the solid will be polarized 
in the plane perpendicular to the plane of incidence. If the incident light is polarized 
parallel to the plane of incidence, the intensity of the reflected light will be theoret- 
ically zero at the Brewster angle. 

The proposed experiment used the polarization by reflection because of its 
simplicity, but other polarization methods can be used too. When a light fascicle 
falls on the M mirror at a Brewster angle, the reflected fascicle is linearly polarized. 
Using a second rotating mirror, M’, the Malus Law can be checked. If mirror M’ is 
rotated around the PP’ axe, the reflected P’S’ fascicle has a variable intensity with 
two minimum and two maximum values. When the second fascicle falls as the first 
fascicle at the Brewster angle, the S’P’ fascicle has a minimum value (Fig. 4.9). 

The experiments with visible light and mirrors were made more than a century 
ago, starting with Brewster and repeated in a lot of laboratories. To date some 
experiments have been done in order to check the phenomena for IR, and the effect 
can be accepted as valid for UV and gamma rays. Related to this kind of experiment, 
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Fig. 4.9 Polarization by 
reflection 





Direction of 
Propogaton 









Electric Field 
Vector 


Magnetic Field 
Vector 


Plane of 
Polarization 


Fig. 4.10 Definition of a polarization vector 


itis worth mentioning a paper by Lambert et al. from 1969: “Rotatable polarizers for 
the near infra-red” [2]. 

Polarization is an important optical property inherent in all laser beams. Brewster 
windows, reflective phase retarders, and absorbing thin-film reflectors use the 
advantage of polarization. On the other hand, it can cause troublesome and some- 
times unpredictable results when ignored. Because virtually all laser sources exhibit 
some degree of polarization, understanding this effect is necessary to specify 
components properly. The following text gives a basic polarization definition and 
presents the polarization types most commonly encountered. 

Light is a transverse EM wave; this means that the electric and magnetic field 
vectors point perpendicular to the direction of wave’s travel (Fig. 4.10). When all the 
electric field vectors for a given wave-train lie in a plane, the wave is said to be plane 
or linearly polarized. The orientation of this plane is the direction of polarization. 


362 4 The Fundamental of Electrodynamics 


Fig. 4.11 Unpolarized light 


Fig. 4.12 Plane-polarized 
wave that points at 45° to 
the axes 





Unpolarized light refers to a wave collection that has an equal distribution of electric 
field orientations for all directions (Fig. 4.11). Although each individual wave-train 
may be linearly polarized, there is no preferred direction of polarization when all the 
waves are averaged together. 

Randomly polarized light is exactly what it says—the light is plane-polarized, but 
the direction is unknown and may vary with time. Random polarization causes 
problems in optical systems because some components are polarization sensitive. 
If the polarization state changes with time, then the components’ transmission, 
reflection, and/or absorption characteristics will also vary with time. 

Polarization is a vector that has both direction and amplitude. Like any vector, it is 
defined in an arbitrary Coordinate system as the sum of orthogonal components. 
Figure 4.12 shows a plane-polarized wave that points at 45° toward the axes of our 
coordinate system. Thus, when described in this coordinate system, it has equal x and 
y components. If we then introduce a phase difference of 90° (or one-quarter 
wavelength) between these components, the result is a wave in which the electric 
field vector has a fixed amplitude but with a direction that varies as we move down 
the wave-train (Fig. 4.13). Such a wave is said to be circularly polarized because the 
tip of the polarization vector traces out a circle as it passes a fixed point. 

If we have two wave-trains with unequal amplitude and with a one-quarter wave 
phase difference, then the result is elliptical polarization. The tip of the polarization 
vector will trace out an ellipse as the wave passes a fixed point. The ratio of the major 
to the minor axis is called the ellipticity ratio of the polarization. Always state the 
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Fig. 4.13 Circularly 
polarized wave 


polarization orientation when ordering optical coatings for use at non-normal inci- 
dences. If you are unsure about how to determine the polarization state of your 
source, contact applications engineers for assistance. 

When light strikes an optical surface, such as a beam splitter, at a 
non-perpendicular angle, the reflection and transmission characteristics depend on 
polarization. In this case, the coordinate system we use is defined by the plane 
containing the input and reflected beams. Light with a polarization vector lying in 
this plane is called p-polarized, and light that is polarized perpendicular to this plane 
is called s-polarized. Any arbitrary state of input polarization can be expressed as a 
vector sum of these s and p components. To understand the significance of s and 
p polarizations, examine the graph that shows the single surface reflectance as a 
function of angle of incidence for these components of light at a wavelength of 
10.6 um striking a ZnSe surface (Fig. 4.14). Note that although the reflectance of the 
s component steadily increases with angle, the p component at first decreases to zero 
at 67° and then increases after that. The angle at which the p reflectance drops to zero 
is called the Brewster’s angle. This effect is exploited in several ways to produce 
polarizing components or uncoated windows that have no transmission loss such as 
the Brewster windows. 

The angle at which p reflectance drops to zero (i.e., Brewster’s angle) can be 
calculated from 63 = tan '(n), where Og is the Brewster’s angle and n is the 
material’s index of refraction. Polarization state is particularly important in laser- 
cutting application as well as a consideration for purposes of this book and laser 
interaction with metal. Figure 4.15 shows a summary of polarization. 

In addition, by now we know that if the waves travel down the path of say the 
x-direction, it is called transverse because the displacement (in the case of rope or 
string) is perpendicular to the direction of propagation. it is known as longitudinal if 
displacement from equilibrium is along the direction of propagation. Sound waves, 
which are nothing but compressed waves in air, are longitudinal although EM 
waves, as we will see, are transverse (see Sect. 6.2 of this book). 
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Fig. 4.14 For s-polarization, the input polarization is perpendicular to the plane (shown in color) 
containing the input and output beams. For p-polarization, the input polarization is parallel to the 
plane (shown in color) containing the input and output beams 


Polarization Summary 
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orientation. oriented radially to orientation vector 


the beam. describes a circle. 


Fig. 4.15 Polarization summary 
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As we discussed in Sect. 3.2 of this book, we now pay attention to sinusoidal waves 
of frequency w, where in visible range, each unlike frequency corresponds to 
different colors. Such waves are called monochromatic waves (see Table 4.2). 
Suppose, furthermore that the waves are traveling in the z-direction and have no 
x or y dependencies; these are called plane waves [3] because the fields are uniform 
over every plane perpendicular to the direction of propagation (Fig. 4.16). 




























































































Table 4.2 The Frequency (Hz) Type Wavelength (m) 
electromagnetic spectrum I 022 10 13 
107! Gamma rays 10° '? 
102° 197!! 
10!° 19~!° 
10'8 X-rays 10°° 
1017 10-8 
10'° Ultraviolet 10-7 
10° Visible 10° 
10'* Infrared 10° 
103 10-4 
10"? 10° 
10!! 10-7 
10° Microwave 10°! 
10° 1 
10° TV, FM 10 
10’ 10° 
10° AM 10° 
10° 104 
104 RF 10° 
10° 10° 
The visible range 
Frequency (Hz) Color Wavelength (m) 
1.0 x 10" Near ultraviolet 3.0 x 10-7 
7.5 x 104 Shortest visible blue 4.0 x 1077 
6.5 x 10!4 Blue 4.6 x 1077 
5.6 x 10'4 Green 54x 10’ 
5.1 x 10!4 Yellow 5.9 x 107-7 
4.9 x 10'4 Orange 6.1 x 10-7 
3.9 x 10'4 Longest visible red 7.6 x 1077 
3.0 x 10!4 Near infrared 1.0 x 10° 
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Fig. 4.16 Plane wave x 
propagation [1] | 





We are interested then in complex form of both magnetic and electric fields as 
follows: 


E(z,t)= Eye and B(z,t) = Bye” (4.22) 


where E 0 and B gare the complex amplitudes [1]. The physical fields, of course, are 


the real parts of E and B. Now, the wave equations for E and B (see Eq. 4.2) were 
derived from Maxwell’s equations. Still, whereas every solution to Maxwell’s 
equations (in empty space) must obey the wave equation, the converse is nor true; 


Maxwell’ S equations i impose extra constraints on E 9 and B o. In particular, because 
V-E=Oand V-B= 0, it follows that: 


(Eo). = (Bo) =0 (4.23) 


This is because the real part of E differs from the imaginary part only in the 
replacement of sine by cosine; if the former obeys Maxwell’s equations, so does the 


latter, and thus E as well [1]. Equation 4.23 also indicates that EM waves are 
transverse—that is, the electric and magnetic fi fields are perpendicular to the direction 
of propagation. Moreover, Faraday’s Law, VxE=—0B / Ot, implies a relation- 
ship between the electric and magnetic amplitudes in the following equation: 


—k(Eo), = @(Bo), and —k(Eo), = @(Bo), (4.24) 
or in compact form complex variable presentation we have: 


(2x Eo) (4.25) 


From this equation it is evident that E and B are in phase and mutually perpendic- 
ular; their (real) amplitudes are related by 
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Fig. 4.17 Scalar product of c k 


k-r » a 





Bo = —Eo = — Fo (4.26) 


The fourth of Maxwell’s equations, Vx B= peo (OE / or), does not yield an 
independent condition; it simply reduces to Eq. 4.24. 

We can easily generalize to monochromatic plane waves traveling in an arbitrary 
direction. The notation is facilitated by the introduction of propagation (or wave) 
vector k pointing in the direction of propagation, with a magnitude that is the wave 
number k. The scalar product, k-#, is the appropriate generalization of kz in 
Fig. 4.17, so we have: 


7 : (i . ies : (4.27) 
B(#,t) =—Epe'\*?-) (EK x #) =-(k XE 
(7,1) = Lyell) (Ex) = 1% x B) 
where 77 is the polarization vector, and because E is transverse, 
ni-k =0 (4.28) 


The transversality of B follows automatically from Eq. 4.27. The actual (real) 
electric and magnetic fields in a monochromatic plane wave with propagation vector 
k and polarization 7 are [1]: 


=> 


E(?,t) = Ey cos (k- F-—ot+6 


wm 
ay 


(4.29) 
B(rt) = * Fy cos (k-#- wt +6)(k x 7) (4.30) 


Example 4.2 If E points in the x-direction, then B points in the y-direction Eq. 4.25. 
See Fig. 4.18. 
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Fig. 4.18 Depiction of 
Example 4.2 solution 








Fig. 4.19 Sketch of the Example 4.3 solution 


Solution Writing the following 
= ae nae A = ls f(z A 
E(z, t) = Eyeil&7-#) z and B(z, t) hy, Eyeil#r-a) y 
Cc 


or taking the real part, we have 
E(z,t) = Ey cos (kz — wt + 5)X 
7m 1 . (4.31) 
B(z,t) =—Eo cos (kz — of + 6) ¥ 

c 


Example 4.3 Write down the (real) electric and magnetic fields for a monochro- 
matic plane wave of amplitude Eo, frequency w, and phase angle zero that is 
(1) traveling in the negative x-direction and polarized in the z-direction; (2) traveling 
in the direction from the origin to the point (1, 1, 1), with polarization parallel to the 
xz-plane. In each case, sketch the wave and give the explicit Cartesian components of 
kK and it. See Fig. 4.19. 
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Solution We can start with the following for part (a): 
@ k=-2%, a=7 k-¥= (-22)-@¥ +yF +22) =-S x 
c c c 
kxm=-¥xZ=¥ 
za @ 4 op Eo @ + 
E(x, t) = Eo cos (2x4 ot) Z B(x,t) =—cos (=: + wt) y 
c c c 


and we have the following for part (b): 


> 


(b) k= 2 (54) i= unas (Because ii is parallel to the xz-plane, it must 


c V3 v2 
have the form ax + BZ; because i-k =0,fB=—a; and because it is a unit 
vector, a = 1/ J/2): 
> Ou, 4 4 * * @ 
kA PS le ey He) (ae yy ez) See) 
q+ +2): (OX + y¥ +22) Fee ty te) 


~p 
x 
=p 
I 





E(x,y,z, t) = Eocos rete t 9 z)— of (: a ‘ 


~ E 4 3 3 
B(x, y, z,t) = alae ates vs z) of ( X+2y ‘ 
c 











V3¢ 


4.7 Boundary Conditions: The Reflection and Refraction 
(Transmission) Dielectric Interface 


The important problem of the reflection and refraction (transmission) of a wave at the 
interface of two media of different dielectric constants is a boundary value one in 
which the physics is identical in principle to that involved in time-harmonic prob- 
lems involving conductors. If, for example, one medium is vacuum, and a plane 
wave is incident on a second (dielectric) medium, the incident harmonic of the plane 
wave generates oscillating time-harmonic dipoles (or dipolar currents) that produce a 
field of their own [4]. 
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The strength of the currents is not known in advance, however, and therein lays the 
essence of the problem. Although this is a boundary value problem, it is especially 
simple because of the great symmetry and simple geometry. It turns out that it can be 
solved by adding to the incident plane waves only two other plane waves—a reflected 
and a transmitted (or refracted) one. The geometry is shown Fig. 4.20 with the plane 
z = 0 taken as the interface between the two media, labeled 1 and 2. 

In terms of the vector propagation constant, Ki, the incident fields are 


E; = Ae’? and B; = (Ki/k) x Ej Incident (4.32) 


Here B; is derived from the assumed EF. iby V x E i= ikB;. The normal to the plane 
z = 0 and the vector k; defines a plane of incidence that, without loss of generality, 
can be taken to be the x, z-plane as was shown in Fig. 4.18. We now postulate the 
existence of two other plane waves and will show that these suffice to solve the 
boundary value problem. These reflected and transmitted waves, amplitudes E, and 
E,, respectively, are [4]: 


m= 


E, = Re’? and B= (k,/k) x E, Reflected (4.33) 
E,—Te*? and B, = (k,/k) x E, Transmitted (4.34) 


The vectors k, and ki, for the moment, must be considered arbitrary in direction, 
for although k; in the x, z-plane, we cannot assume the same a priori for k, and k;. The 


magnitudes, k; = 





k and so on, of the wave vectors are, with 1; and nz the indices of 


refraction of the two media and k = a/c, 


kj =k, = nk 


4.35 
k; => nok ( ) 


Considering the boundary condition that states the tangential components of E 
and B are continuous, we have to satisfy the following two steps: 


First: If the tangential components of the three fields in Eqs. 4.33, 4.34, and 4.35 are 
to be matched at z = 0, it also should be clear that the spatial dependence given by 
the exponents must be identical. However, this is a necessary but not a sufficient 
condition. 


Second: Vector coefficients ZA: R, and T must be determined. 


The first condition—that is, the spatial variation of the three fields must be identical 
at z = 0—leads to the following equation: 


(ki , tax = (k, ' Fy 2%, = (k, , es (4.36) 


The first equality in Eq. 4.35 yields kx = k,,.x + k,,y. For this condition to hold 
for all x and y, we must have k,,, = 0, showing that k, lies in the plane of incidence 
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and, also that k;,.x = KX. Similarly, from the second equality in Eq. 4.36, k, must lie 
in this plane so that i, ky, and k; are all co-planar. Moreover, from the geometry of 
Fig. 4.19, we have from k;, = k,, that k; sin 0; = k, sin @,. Or, because of k; = k,, we 
can write the following: 


sin6; = sin 6, (4.37) 


The angle of incidence equals the angle of reflection. Similarly, the equality of k;x 
and k,, yields: 





k; sin 0; = k, sin 0; (4.38) 
Or, using Eq. 4.35, we have: 
sin; no 
= 4.39 
sin 0; ny ( ) 


Equation 4.38 is known as Snell’s Law of refraction. The condition of Eqs. 4.37 
and 4.38 are quite general ones that are independent of the detailed vectorial nature 
of the wave field. They hold for reflection and refraction of scalar waves (SWs). 
These conditions by themselves do not guarantee the continuity of tangential E and 
B across the boundary. To satisfy these conditions more specifications of the 
polarization of the fields should be analyzed [4]. For convenience we consider the 
general case of arbitrary incident polarization as a linear combination of a wave with 
polarization perpendicular to the plane of incidence and one with polarization 
parallel to this plane. The reflected and transmitted waves will then be similarly 
polarized. These two cases are sketched as parts (a) and (b) in Fig. 4.20. 

Consider first the case of E perpendicular to the plane of incidence (i.e., in the 
y-direction). Vector coefficients A, R, and T become scalar ones, with a subscript L 
to denote this case: 





Fig. 4.20 Reflection and refraction at a dielectric interface: (a) Electric vector perpendicular to 
plane of incidence. (b) Electric vector parallel to plane of incidence 
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Ey, = eke 7 
E,, = elke r (4.40) 
E, = i k, r 
Then from the continuity of tangential E at the boundary, we have 
AL +R, =T. (4.41) 
The condition on tangential B becomes, with B;, = —A, cos 0;, B,. = —R, cos 0,, 
and B, = —T, cos 0,—also considering that cos?; = —cos 0,—we then have: 


n(A, — R,)cos6; = mT , cos0, (4.42) 


Solving Eqs. 4.41 and 4.42 for the ratios R,/A, and T/A, we find the following 
using Snell’s Law: 

















tan 0; 
R, tan, sin (0; — 0;) 
i 1 4 andi ~~ sin (6; + 8) 
tan 0, (4.43) 
T 2 2 cos 9; sin 0, 
Ay, , tan ~ sin (6; + 4) 





tan 0, 


For the second case, an incident wave polarized parallel to the plane of incidence, 
we use A) = | A| in Eq. 4.33 and, similarly, R) = [R| and 7 = | T| in Eqs. 4.34 and 
4.35. The boundary condition yields: 


cos 0; (Aj _ Ri) = cos 6,7 | 








(4.44) 
ny (Aj — Ry) = mT 
These equations lead to the results: 

R\ _ tan (0; — 4) 
A tan (0; + 6) 

: (4.45) 
T| 7 2 cos 0; sin 0, 
Aj sin (6; + @,) cos (8; — @,) 





There are two phenomena worthy of note in connection with the preceding 
discussion. Consider the case of polarization in the plane of incidence. We can see 
from Eq. 4.43 that R\/A) will be zero for 0; + 0, = 2/2. Putting this condition into 
Snell’s Law, using sind, = sin (2/2 — 6;) = cos 6;, we can see that the angle of 
incidence, Og (i.e., the Brewster angle) for which this happens, is defined by 


4.7 Boundary Conditions: The Reflection and Refraction (Transmission). . . 373 


tan 03 = (4.46) 
ny 

If a wave with arbitrary polarization is incident on a dielectric interface, it can be 
considered to be a linear combination of a wave polarized parallel to, and a wave 
polarized perpendicular to, the plane of incident. At the Brewster’s angle, the parallel 
component will not be reflected, so the reflected wave will be plane-polarized in a 
plane perpendicular to the plane of incidence. This effect can then be made the basis 
of a device for polarizing an unpolarized beam of radiation. 

The second phenomenon is that of total internal reflection. In either (a) or (b), 
shown earlier in Fig. 4.4, suppose that the index nj, is greater than nz. Then, from 
Snell’s Law, sin@, = (1,/n2) sin 0; is always greater than 6;. There then will be some 
value of 6;, call it Aint, for which 0, = z/2; this angle is defined by 


sin Bint = (4.47) 
n\ 
Because, in general, 
E, = Teiknalx sin 0;+z cos @;) (4.48) 


for 0, = 2/2 there is no wave in the second medium; the z-dependence vanishes. Now 
if 0; > Oint, then sind, is larger than unity from Snell’s Law and, as a consequence, 
cos@, is imaginary, as follows: 


2. 2 
cos 0, = i t= (=) sin20; = 7 (=) sin20; — 1 (4.49) 
n2 n2 


Equation 4.36 then becomes 


E, = Te74/ (n; /ny) sin20;- el sin (4.50) 


This corresponds to a wave that is exponentially attenuated as a function of z, and 
that propagates as a function of x with a propagation constant of kn2 sin 8, Such a 
wave is the prototype of a surface wave; readers can see further discussion and more 
details in Eyges (1972) [4] in his Appendix 2. 

The interesting question is this: What happens when a wave passes from one 
transparent medium into another—air to water, say, or glass to plastic? As in the case 
of waves on a string, we expect to get a reflected wave and a transmitted wave. The 
details depend on the exact nature of the electrodynamics boundary conditions [1]. 





eEt =eEt (i) E!—El (iii) 
1 1 (4.51) 

Bi=By (ii) —B] =—8) Wy) 

My H2 
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These equations relate to the electric and magnetic fields just to the left and just to the 
right of the interface between two linear media. In the following sections we use 
them to deduce the laws governing reflection and refraction of EM waves. 


4.8 Electromagnetic Waves in Matter 


Some solutions to Maxwell’s equations already have been discussed in previous 
sections. This section extends the treatment of EM waves. Because most regions of 
interest are free of charge, it will be assumed that charge density p = 0. Moreover, 
linear isotropic (invariant with respect to direction) materials will be assumed, with 
the following relationships: 


D=cE 
B=yH (4.52) 
J=cE 


where 


D = flux density (C/m?) 

E = electric field (N/C) 

€ = permittivity of the medium (C?/N ‘ m?) or equivalent (F/m) 
B= magnetic field (T) 

H= magnetic field strength (A/m) 

ft = mobility within materials (m7/V - s) 

J = current density (A/m?) 

o = conductivity of materials (S/m) 

p = charge density (C/m*) 


where 

A = Ampére 
C = Coulomb 
N = Newton 
F = Faraday 
S = Siemens 
T = Tesla 


4.8.1 Propagation in Linear Media 


Inside matter, but in regions where there is no free charge or free current, Maxwell’s 
equations become 
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4 5 B 
(i) V-D=0 (iii) vxi=-2 
5B (4.53) 
- " D 
(ii) V-B=0 (iv) Vita 
If the medium is Jinear, then 
= 2 fs 1s 
D=eE and H = a (4.54) 


and homogeneous (so ¢ and yz do not vary from point to point), Maxwell’s equations 
reduce to: 


(i) V-E=0 (iii) vxE=- 
7 (4.55) 
(ii) V-B=0 (iv) Vx Ba pee 


which remarkably differ from the vacuum analogs in Eq. 4.1 only in the replacement 
of Hoéo by pe. This is obvious mathematically, yet the physical implications are 
astonishing [5]. As the wave passes through the fields, they busily polarize and 
magnetize all the molecules, and the resulting (oscillating) dipoles create their own 
electric and magnetic fields. These combine with the original fields in such a way as 
to create a single wave with the same frequency but a different speed. This extraor- 
dinary conspiracy is responsible for the phenomenon of transparency. It is a 
distinctly non-trivial consequence of the linearity of the medium [5]. Looking at 
the set of Eq. 4.54, it is evident that EM waves propagate through a linear homoge- 
neous medium at a speed: 


1 Cc 
= = 4.56 
ae (4.56) 


ee (4.57) 
€0Ho 


is the index of refraction of the material. For most materials yz is very close to fg, so 


n& fe, (4.58) 


where ¢€, is the dielectric constant, also known as relative permittivity [1], and is 
equal to €, = e/e9. Because e, almost always is greater than 1, light travels more 
slowly through matter—a fact that is well known from optics [1]. 


where 


n 
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Fig. 4.21 Normal wave 3 
incidence [1] 
E; © 


v2 
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4.8.2 Reflection and Transmission at Normal Incidence 


Suppose the xy-plane forms the boundary between two linear media. A plane wave 
of frequency @, traveling in the z-direction and polarized in the x-direction, 
approaches the interface from the left (see Fig. 4.21): 


E,(z, t) = Eo, e200) = 


B 1 ey pilkiz—on) 2 (4.59) 
By(z,t) = Fy eee ¥ 
vi 
This gives rise to a reflected wave 
Ex(z, th = Ep, eis X 
B LS Meas (4.60) 
Br(z, t) = — Fi, ¥ 
v1 
which travels back to the left in medium (1), and a transmitted wave 
Er(z, t) — Eo, elthiz- 9 3 
(4.61) 


3 Le. saree 
Br(z, t) = — Boe y 
1 


which continues to the right in medium (2). Note the minus sign in Bp, as required by 
Eq. 4.27, or, if you prefer, by the fact that Poynting’s vector aims in the direction of 
propagation. 7 7 7 7 

At z = 0 the combined fields on the left, E; + Ep and B, + Br, must join the 


fields on the right, Er and By, in accordance with the boundary conditions in 
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Eq. 4.50. In this case there are no components perpendicular to the surface, so (i) and 
(ii) are trivial. However, (iii) requires: 





Eo, + Eo, = Eo, (4.62) 
while (iv) says 
I fle 1 ~ 1/l- 
nar ee ee 8) 
Hy \Y1 v\ Hz \02 
or 
Eo, — Eo, = BEo, (4.64) 
where 
— H1%1 _ Hila (4 65) 
#202 M2M1 


Equations 4.61 and 4.63 are solved easily for the outgoing amplitudes in terms of the 
incident amplitude: 


e 1- ws ws 2 ~ 
op = (5) and Eo, = (5) (4.66) 


These results are strikingly similar to the ones for waves on a string. Indeed, if the 
perrnittivities 4 are close to their values in vacuum (as, remember, they are for most 
media), then 6 = *, and we have [1]: 


~ D2 — V1 \ x ~ 2v2 zs 
Eo, = E d £o, = E 4.67 
‘ (2 =") Ga Be (; "| ‘ (4.67) 


In that case, as before, the reflected wave is in phase (right side up) if v2 > v; and out 
of phase (upside down) if v2 < v1; the real amplitudes are related by 














2s 1—p 2v2 
Or (; ze ,) an Or (— Eo -) 0; ( ) 
or, in terms of the indices of refraction [1], 
E mT \ eo, and E, 2m 5 (4.69) 
= an = | ——— : 
Oe ny + Ng " = ny +n Mm 


According to the definition of intensity (average power per unit area) that is given 
by Griffiths [1], as follows, we should ask these two questions: Which fraction of the 
incident energy is reflected? Which fraction is transmitted? 


1 
I= 580k (4.70) 
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If again 4 = #2 = Uo, then the ratio of the reflected intensity to the incident intensity 


is as follows: 
1 E,.\2 = 2 
a ( 2) = (= =) (4.71) 
Tr Eo, ny +n 


whereas the ratio of the transmitted intensity to the incident intensity is: 





T= 





Iv _ E202 2) _ 4nyn2 
( 


4.72 
Eo, ny t+ m) ( ) 


Ty eyvy 


R is called the reflection coefficient and T is called the transmission coefficient; 
they measure the fraction of the incident energy that is reflected and transmitted, 
respectively. Notice that [1]: 


R+T=1 (4.73) 


as conservation of energy, of course, requires. For instance, when light passes from 
air (n; = 1) into glass (nz = 1.5), then R = 0.04 and T = 0.96. Not surprisingly, most 
of the light is transmitted [1]. 


Example 4.4 Calculate the exact reflection and transmission coefficients, without 
assuming fy = 2 = Mo. Confirm R + T= 1. 
E R 


2 
Solution From Eq. 4.65 we have R = (=) ; substituting from Eq. 4.65 results, 
I 


2 
where f = ae and using Eq. 4.66 results in T = (35) , which is Eq. 4.60. Note 


the following: 





2 
£202 fly E202 02 fly (>) Wy 
€1V) Ho €101 vi Ho \U2 


202 
R= +p: [46+ (1 - 6) | 
=a get 1 — 26+ p’) 


= appl tte =1 








Example 4.5 In writing Eqs. 4.60 and 4.61, the assumption was that the 
reflected and transmitted waves have the same polarization as the incident 
waves—along the x-direction. Prove that this must be so. [Hint: Let the polarization 
vector of the transmitted and reflected waves be 7i = cos Ore + sin Ory and i R= 
cos rx + sin Ory, and prove from the boundary condition that 07 = 0z = 0.] 
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Solution Equation 4.52 is replaced by Eo,x + Eo, it a aa and Eq. 4.58 by 
Eo,¥ - Eo, (Z x rir) = Eo, (Z x ir). The y-component of the first equation is 
E Op SINOR = E 0, Sin @7; the x-component of the second is E Op SINOR = —fBEo, sin Or 
Comparing these two, we can conclude that sindg = sin 067 = 0, thus, 07 = Or = 0. 


4.8.3 Reflection and Transmission at Oblique Incidence 


In the previous section reflection and transmission were treated at normal inci- 
dence—that is, when the incoming wave hits the interface head-on at 90° incident 
angle [1]. The more general case of oblique incidence, in which the incoming wave 
meets the boundary at an arbitrary angle 6), also is treated by David Griffiths (see 
reference [1] of this chapter). We just wrote the result of his conclusion here and 
encourage readers to refer to his book’s Sect. 9.3.3 (p. 386) for more details. 

Although the angle of incidence 0; = 0 is a special case of the oblique incidence, 
we did treat that separately for some cases; here we are dealing with high-power laser 
interactions with matter. The treatment by Griffiths is based on a monochromatic 
plane wave with an arbitrary angle of incidence (Fig. 4.22). 

Suppose that a monochromatic plane wave (incident wave) approaches from the 
left as follows: 


a > (£7 aie 1 > = 
E,(7,t) = Bo,e! (tr For) and B, (7 t) =—(k x E;) (4.74) 
This gives rise to a reflected wave of the form in Eq. 4.69: 


E,(?, t) = Eo, ei (fe ?-o1) and Br (7, t) =— (Kr x Ep) (4.75) 


Fig. 4.22 Incidence at 
oblique angle [1] kr 





N 






Plane of Incidence 
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and a transmitted set of wave equations of: 


Er (7,1) = Eo,e() and By (Ft) =—(kr x Er) (4.76) 
All three waves have the same frequency, @, that is determined once and for all at the 
source (e.g., laser beam, flashlight, or whatever) that produces the incident beam. 
The three wave numbers are related per Eq. 4.11: 

v2 


ko bo ow oF hee ee oe (4.77) 
vi nN 


The combined fields in medium (1), E r+ E rand B r+ Br, must now be joined to 


the fields E'; and By in medium (2), using the boundary conditions in Eq. 4.51. These 
all share the generic structure. Boundary conditions should hold at all points on the 
plane, and for all times for all three plane-wave equations when z = 0. Otherwise, a 
slight change in x would destroy the equality (see Example 4.6) [1]. The time factors 
in this case already are equal; in fact we can regard this as an independent confir- 
mation that the transmitted and reflected frequencies must match the incidence 
frequency. So, for a spatial case, we have the following relation: 


ky ?=kp-?=kp-? when z=0 (4.78) 
or more explicitly, 
x(kr), + y(kr)y = x(kr)» + y(kr)y = x(kr), + y(kr)y (4.79) 


for all x and y. Equation 4.73 can only hold if the components are separately equal, 
for if x = 0, we get: 


(kr), = (kr), = (kr), (4.80) 
While y = 0 gives 


(ki)x = (kr)x = (kr) x (4.81) 


We may as well orient the axes so that k lies in the xy-plane (i.e., (k7) = 0); according 
to Eq. 4.74, so too will kr and kr. 


Conclusion 

First Law: The incident, reflected, and transmitted wave vectors form a plane (called 
the plane of incidence), which also includes the normal to the surface (here, the 
z-axis). Meanwhile, Eq. 4.81 implies the following: 


ky sin 0; = kr sin Or = kr sin Or (4.82) 


where 6; is the angle of incidence; Op is the angle of reflection; and @, is the angle of 
transmission, more commonly known as the angle of refraction—all of them are 
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measured with respect to normal (see Fig. 4.21). In view of Eq. 4.81, we get the 
second law. 
Second Law: The angle of incidence is equal to the angle of reflection, 


0, = Or (4.83) 


This is the law of reflection. 
Third Law: The transmitted angle follows: 


sin Or ny 





sin@; nm Ge) 
This is the law of refraction, or Snell’s Law. 

These are the three fundamental laws of geometrical optics. It is remarkable how 
little actual electrodynamics went into them; we have yet to invoke any specific 
boundary conditions. Now that we have taken care of the exponential factors they 
cancel, given Eq. 4.77 the boundary conditions in Eq. 4.51 become [1]: 





(i) €1(Bo, + Eo), = €2(Eor) 

(ii) (Bo, + Bo,), = (Bor), 

(ii) (Ba, + Bx), = (Eo), ay 
(iv) 7(By, + Boy),, = (Bo), 


where By = (1/v) k x Eo in each case. The last two, (iii) and (iv), represent pairs of 
equations—one for the x-direction and one for the y-direction. 


Suppose that the polarization of the incident wave is parallel to the plane of 


incidence (the xz-plane in Fig. 4.23); it follows that the reflected and transmitted 


Fig. 4.23 Incident 
oblique wave 
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waves also are polarized in this plane (see Example 4.5 in this chapter). We present 
this in an sample to analyze the case of polarization perpendicular to the plane of 
incidence (see Example 4.7). Then (i) reads: 


€\(—Eo, sin 0; + Eo, sin Op) = €r(— Eo, sin @r) (4.86) 


Note that (ii) adds nothing (0 = 0) because the magnetic fields have no z-compo- 
nents; (iii) becomes [1]: 


Es, cos 6; + Ey, COS Or = Eo, cos Or (4.87) 
and (iv) says, 
Nie a el ae (4.88) 
HyV1 M202 


Given the laws of reflection and refraction, Eqs. 4.84 and 4.88 both reduce to the 
following: 





Eo, — Eon = BEo, (4.89) 
where again 
— #11 _ Mile (4 90) 
#202 MQM 
Equation 4.87 says: 
Eo, + Eo, = Eo, (4.91) 
where again 
cos Or 
= 4.92 
” cos 0; ( ) 


Solving Eqs. 4.89 and 4.91 for the reflected and transmitted amplitudes, we 
obtain [1]: 





(4.93) 





a 2 a 
Eo, = (=) Eo, 
a+p 


These are known as Fresnel’s equations, for the case of polarization in the plane of 
incidence. (There are two other Fresnel equations that give the reflected and trans- 
mitted amplitudes when the polarization is perpendicular to the plane of incidence— 
see Prob. 9.16 in Griffiths [1].) Notice that the transmitted wave is always in phase 
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with the incident one; the reflected wave is either in phase (right side up) if a > f, or 
180° out of phase (upside down) if a < f. 

Be aware that there is an unavoidable ambiguity in the phase of the reflected 
wave. The convention that Griffiths [1] used in Chap. 9 of his book, which is 


reflected here and was adopted for Fig. 4.23, has Ep as positive “upward”; this is 
consistent with some, but not all, of the standard optics’ texts. Changing the sign of 
the polarization vector is equivalent to a 180° phase shift. 

The amplitudes of the transmitted and reflected waves depend on the angle of 
incidence, a, because a is a function of 0;: 











L— sin26r _ i 1 — [(m1/np) sin 0,)” (4.94) 


cos 6; cos 0; 


In the case of normal incidence (0; = 0), a = 1 so we retrieve Eq. 4.60. At grazing 
incidence (0; = 90), a diverges, and the wave is totally reflected (a fact that is 
painfully familiar to anyone who has driven at night on a wet road). Interestingly, 
there is an intermediate angle, O03 (i.e., Brewster ’s angle), at which the reflected wave 
is completely extinguished [1]. 

Because waves polarized perpendicular to the plane of incidence exhibit no 
corresponding quenching of the reflected component, an arbitrary beam incident at 
Brewster’s angle yields a reflected beam that is totally polarized parallel to the 
interface. That is why Polaroid glasses, with the transmission axis vertical, help to 
reduce glare off a horizontal surface [1]. 

According to Eq. 4.94, this occurs when a = f, or 


1-# 
sin?@g = —— (4.95) 
(n/m) — B 
For the typical case py, & po, so B & no/nj, sin’?Og & B7/(1 + fH”); thus, 
tan 0g & ue (4.96) 


Figure 4.24 shows a plot of the transmitted and reflected amplitudes as functions 
of 0;, for light incident on glass, (n2 = 1.5), from air, (nj = 1.0). (On the graph, a 
negative number indicates that the wave is 180° out of phase with the incident 
beam—the amplitude itself is the absolute value). 

The power per unit area striking the interface is S- 2. Thus, the incident 
intensity is: 


1 
I; = 5 £1r1£5, cos Or (4.97) 


while the reflected and transmitted intensities are 
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Fig. 4.24 Plot of the transmitted and reflected amplitude [1] 
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Fig. 4.25 Plot of reflection and transmission coefficient for polarized waves 


1 2 1 2 
I; = 51M Ep, cos@; and Ir = 5 e2v2Eo, cos Or (4.98) 
(The cosines are there because we are talking about the average power per unit area 
of interface, and the interface is at an angle to the wave front [1].) 

The reflection and transmission coefficients for waves polarized parallel to the 
plane of incidence are 








ex)’ (: ; Al 
1 (= a +f ( ) 
Ir €2U2 Eo 2 COs Or 2 2 
fos = z = ap{ —— (4.100) 
I; E,V1 Eo, cos 0; a+ fp 


They are plotted as functions of the angle of incidence in Fig. 4.25 (for the air—glass 
interface). R is the fraction of the incident energy that is reflected—naturally, it goes 
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to zero at Brewster’s angle; and T is the fraction transmitted—it goes to | at Og. Note 
that R + T = 1 is required by conservation of energy; the energy per unit time 
reaching a certain patch of area on the surface is equal to the energy per unit time 
leaving the patch [1]. 

1ax 


ibx __ 


Example 4.6 Suppose Ae“ + Be’”* = Ce‘, for some non-zero constant A, B, C, a, 
b, c, and for all x. Prove thata = b=candA+B=C. 


Solution Ae“ + Be’”* = Ce‘, for all x, so (using x = 0), A + B = C. Differentiate: 
iaAe + ibBe’”* = icCe', so (using x = 0), aA + bB = cC. Differentiating again: 
—ad’ Ae — b*Be'”* = —c*Ce'™, so (using x = 0), @A + b°B= °C. a°A+b°B=c 
(cC) = c(aA + bB), or (A + B)(a’A + b?B) = (A + B)c(cC) = c(aA + bB) = cC 
(aA + bB), or a°A + b°AB + a°AB + b’B? = (aA + bB)Y = aA” + 2abAB + b’B’, or 
(a? + b? — 2ab)AB = 0, or (a — b)’AB = 0. 

But A and B are non-zero, so a = b. Therefore, (A + Be = Ce’ - a(A + B) =cC, 
or aC = cC, so (because C # 0), then a = c. Conclusion: a = b = c is driven. 


Example 4.7 Analyze the case of polarization perpendicular to the plane of 
incidence (i.e., electric fields in the y-direction (see Fig. 4.22). Impose the boundary 
condition in Eq. 4.85, and obtain the Fresnel equations for Eg and Eo. Sketch 
(Eo,/Eo,) and (Eo, /Eo,) as functions of / = n/n, = 1.5. (Note that for # the 
reflected wave is always 180° out of phase.) Show that there is no Brewster’s angle 


for any n, and no: Eo is never zero (unless, of course, 1, = nz and py = M2, in which 
case the two media are optically indistinguishable). Confirm that your Fresnel 
equations reduce to the proper forms at normal incidence. Compute the reflection 
and transmission coefficients, and check that they add up to 1. 


Solution We start with the following relationships: 
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(i) Et =e,Et (iii) E! = E! 
Boundary conditions: 1 1 
(i) BE =By (iv) — Bl =a 
faa My 
sin 0 v2 


Law of refractions: 





sin; ov, 


(Note: ky -?-ot=ke-F—ot= ie -F—ot at z = 0, so we can drop all 
exponential factors in applying the boundary conditions.) 


Boundary condition (i): 0 = 0 (trivial) 
Boundary condition (iii): Eo, + Ep = Eo, 
Boundary condition (ii): 1 Eo, sin 0, + 1 Eon sin0, = 1 £p, sin0, = Eo, + Eo, = 


v, sin@> \ F 
(2 sin a) Eo, 
But the term in parentheses is 1, by the law of refraction, so this is the same as (ii). 


Boundary condition (iv): me [2 Eo,(—cos A:) + ao (—cos 6)| = ap Er (—cos 0) 








Hav? sin A, cos 0; 
have Eo, _ Eo, = apEo,. Solving for Eb. and Ey. 2Eo, = (1+ ap)Eo,+ or 
Eo, — (33) Bo, and Eo, = Eo, — Eo, — ( 2 1226) Boy; then we have 


I+ap 1+ap 
Eq, = [— 


Because a and f are positive, it follows that 2/(1 + af) is positive; thus, the 
transmitted wave is in phase with the incident wave, and the (real) amplitudes are 


related by Ep, = (23) E,. The reflected wave is in phase if af < 1 and 180° out of 


oe = sin 02) = s = 
or we can write Eo, — Eo, = (2 = 2) Eo,. Let a = 2°% and p = foe Then we 
‘202 








SS 
SS 
py 

— 


1+ap 
Fresnel equations for polarization perpendicular to the plane of incidence. 


phase if af < 1; the (real) amplitudes are related by Eo, = (53) Eo,. These are the 
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To construct the graphs, note that af = 6 ¥225=s"*@ where @ is the angle of 


cos@ 


vi—sin?d _ V f—sin20 


cos@ cos? 





incidence, so for # = 1.5 and af = fp 











0° 10° 20° 30° 40° 50° 60° 70° 80° 90° 6; 


Is there a Brewster’s angle? Well, Zp, = 0 would mean that af = 1, and thus that 


psec Le ee a (2) sin2o = (222) cos so l= (2) 
Dv, 2 


= cos 6 ~ BT poy? Hy v| 
[sin70 + (u2/p1) cos 76]. 

Because pl, © fo, this means that 1 = (v2/v1), which is only true for optically 
indistinguishable media in which case there is, of course, no reflection—but that 
would be true at any angle, not just at a special “Brewster’s angle.” If #2 were 
substantially different from j;, and the relative velocities were just right, it would be 


possible to get a Brewster’s angle for this case, at 











= cose = (o1/o2)" — (uae2/my61) —1 
(Hi /H2) — 1 (Ha/Hy) — 1 
_ (€2/€1) = (2/m1) 
(Ha/H1) — (Hr /H2) 


But the media would be very peculiar. By the same token, 6p is either always 0, or 
always z, for a given interface—it does not switch over as you change @ the way it 
does for polarization in the plane of incidence. In particular, if # = 3/2, then af > 1, 


for af = (“s) > 1 if 2.25 — sin?0 > cos?0, or 2.25 > sin20 + cos’@ = 1.0. 


cos@ 
In general, for # > 1, af > 1, and thus 6g = a for f < 1, af < 1, and dr = 0. 
At normal incidence, a = 1, so Fresnel’s equations reduce to Eo, = 35) Eo, and 


Eqn = (5) Eo, consistent with Eq. 4.68. 
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Reflection and transmission coefficients: R = (=) 
1 





2 
Referring to Eq. 4.100, T = 28a(z) = ap (35) , therefore R+T= 


e10| Eo, 





(l-af)’+4af _ 1-2ap+a?p?+4ap _ (It+ap)* _ 1 
(tap) (4 ap)” (Lp) 








Example 4.8 The index of refraction of a diamond is 2.42. Construct a graph 
analogous to Fig. 4.23 for the air—-diamond interface (assume fy = 2 = Mo). In 
particular, calculate: (1) the amplitudes at normal incidence, (2) Brewster’s angle, 
and (3) the “crossover” angle at which the reflected and transmitted amplitudes are 
equal. 


“* 2 
Solution In Eq. 4.84 we see that 6 = 2.42 and in Eq. 4.88 a = W-Gnehey 


cos@ 


@6=0 5 a = 1. See Eq. 4.97 (=) anf — 1-242 _ 12 — 0,415, 











(=) = 2, = 2, = 0.585. For Eq. 4.90 +Ogtan~'(2.42) = 67.5° 

Eo, = Eo, > @—B = 2;a=fh+2=4.42 and (4.42)*cos"6 = 1 — sin*0/(2.42)’; 
(4.42)°(1 — sin?@) = (4.42)? — (4.42)°sin?@ = 1 — 0.171sin70; 19.5 — 1 =(19.5 — 017) 
sin’0; 18.5 = 19.3sin70; sin?@ = 18.5/19.3 = 0.959 > sin 0 = 0.979; 0 = 78.3. 





























4.9 Absorption and Dispersion 


Expressions for the stored energy, energy flow, and power dissipated are derived for 
EM waves in terms of complex permittivity € and permeability for a frequency- 
dispersive absorbing medium. This is shown to be possible when é and jf are known 
functions not only of frequency but also of all the loss factors (e.g., collision 
frequencies, etc.) The derivation is not restricted to media with small losses. 


4.9 Absorption and Dispersion 389 


In physics and electrical engineering, dispersion most often refers to frequency- 
dependent effects in wave propagation. Note, however, that there are several other 
uses of the word “dispersion” in the physical sciences. In the presence of dispersion, 
wave velocity is no longer uniquely defined, giving rise to the distinction of phase 
velocity and group velocity. A well-known effect of phase velocity dispersion is the 
color dependence of light refraction that can be observed in prisms and rainbows. 

Dispersion relations describe the interrelations of wave properties (e.g., wave- 
length, frequency, velocities, refraction index, and attenuation coefficient). Besides 
geometry- and material-dependent dispersion relations, there are the overarching 
Kramers-Kronig relations that connect the frequency dependences of propagation 
and attenuation. Dispersion may be caused either by geometric boundary conditions 
(e.g., waveguides, shallow water) or by interaction of the waves with the transmit- 
ting medium. 


4.9.1 Electromagnetic Waves in Conductors 


Previous sections stipulated that free charge density py and free current density Ki if are 
zero, and everything that followed was predicated on that assumption. Such a 
restriction is perfectly reasonable when you are talking about wave propagation 
through a vacuum or through insulating materials such as glass or (pure) water. But 
in the case of conductors we do not independently control the flow of charge and, in 
general, J, is certainly not zero. In fact, according to Ohm’s Law, the (free) current 
density in a conductor is proportional to the electric field: 


Jy =0E (4.101) 


With this, Maxwell’s equations for linear media assume the form 


— - = 8 
(i) V-E=-p, (ili) Vx B= rs 
- (4.102) 
a ae 3 0 
(ii) V- B=0 (iv) V x B= wok + pea 
Now the continuity equation for free charge becomes, 
Sas Op 
Siete 4.103 
f Or ( ) 
Together with Ohm’s Law and Gauss’s Law (i), this gives 
Op > = Oo 
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for a homogeneous linear medium, from which it follows that 
p(t) =e "9 ,(0) (4.105) 


Thus, any initial free charge density pO) dissipates in a characteristic time t = e/o. 

This reflects the familiar fact that if you put some free charge on a conductor, it 
will flow out to the edges. Time constant t affords a measure of how “good” a 
conductor is: For a “perfect” conductor o = oo and t = 0; for a “good” conductor, t 
is much less than the other relevant times in the problem. In oscillatory systems this 
that means t < 1/q; for a “poor” conductor, t is greater than the characteristic times 
in the problem (t > 1/q@) [6]. Ashby points out that for good conductors, it is 
absurdly short (10~!° s for copper, whereas the time between collisions is 7, = 107 '* 
s) [6]. The problem is that Ohm’s Law itself breaks down on time scales shorter than 
T-; actually, the time it takes free charges to dissipate in a good conductor is on the 
order of z,, not tz. Moreover, Ohanian shows that it takes even longer for the fields 
and currents to equilibrate [7]. Nevertheless, none of this is relevant to our present 
purposes; the free charge density in a conductor eventually does dissipate, and 
exactly how long the process takes is beside the point. 

At present we are not interested in this transient behavior—we will wait for any 
accumulated free charge to disappear. From then on p¢= 0, we have 


(i) V-E=0 (iii) vxE=-2 


(4.106) 


a “ xs OE 

(ii) V-B=0 (iv) Vx B= yok + pe~- 

These differ from the corresponding equations for non-conducting media (see 

Eq. 4.55) only in the addition of the last term in (iv). Applying the curl to (iii) and 
(iv), as before, we obtain modified wave equations for E and B: 


3 OE OE OB OB 








VE =e 57 + po om and VB=wpe Ye + yo > (4.107) 
These equations still admit plane-wave solutions, 
B(z,t) = Eye") and B(z,t) = Boel’) (4.108) 
but this time “wave number” k is complex: 
kK? = pew + inow (4.109) 


You can easily check this by plugging Eq. 4.108 into Eq. 4.107. Taking the square 
root, 


k=k+ix (4.110) 
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where 








(4.111) 


The imaginary part of k results in an attenuation of the wave (i.e., decreasing 
amplitude with increasing z): 


E(z, t) = Eoe~Keilke-at) and Bz, th= Boe-Keilke—on) (4.112) 


The distance it takes to reduce the amplitude by a factor of 1/e (about one-third) is 
called the skin depth: 


1 
d== (4.113) 


It is a measure of how far the wave penetrates into the conductor [1]. Meanwhile, the 
real part of k determines the wavelength, the propagation speed, and the index of 
refraction in the usual way: 








l= a, ae (4.114) 


The attenuated-plane waves (see Eq. 4.112) satisfy the modified wave (see 


Eq. 4.107) for any Eo and Bo. Nonetheless, Maxwell’s equations (see Eq. 4.100) 
impose further constraints, which serve to determine the relative amplitudes, phases, 
and polarizations of E and B. As before, (i) and (ii) rule out any z-components; the 
fields are transverse. We may as well orient our axes so that Eis polarized along the 
x-direction: 


E(z,1) = Ee Kell (4.115) 
Then (iii) gives 
B(z, 1) = Ege! leat) § (4.116) 


Part (iv) in Eq. 4.102 says the same thing. Once again, the electric and magnetic 
fields are mutually perpendicular. Like any complex number, k can be expressed in 
terms of its modulus and phase [1]: 


k = Ke'? (4.117) 


where 
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K=|k|=VP +2 =a] eu 1+ (2) (4.118) 


and ¢ = tan '(k/k) 7 - 
According to Eqs. 4.108 and 4.109, the complex amplitudes Ey = Eve’* and Bo 
= Boe'** are related by 
Ke? . 
Boe = ~~ el (4.119) 
@ 





Evidently the electric and magnetic fields are no longer in phase; in fact, 
op — Or = (4.120) 


The magnetic field lags behind the electric filed [1]. Meanwhile, the (real) 
amplitudes of F and B are related by 


al Pate 1+ (2) (4.121) 


The (real) electric and magnetic fields are, finally: 


E(z,t = Eye~™ cos (kz — ot + 5g) X 
(2,1) = Eoe"* | 2) | ne 


B(z, t) = Boe cos (kz — wt + 6g + oy 





These fields are shown in Fig. 4.26. 
Example 4.9 


(a) Suppose you imbedded some free charge in a piece of glass. About how long 
would it take for the charge to flow to the surface? 

(b) Silver is an excellent conductor, but it is expensive. Suppose you were designing 
a microwave experiment to operate at a frequency of 1010 Hz. How thick would 
you make the silver coatings? 

(c) Find the wavelength and propagation speed in copper for radio waves at 1 MHz. 
Compare the corresponding values in air (or vacuum). 


Fig. 4.26 Presentation of x, 
electric and magnetic fields 
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Table 4.3 Resistivities in ohms (All Values Are for 1 atm, 200 °C) [8] 

Material Resistivity Material Resistivity 
Conductors: Semiconductors: 

Silver 1.59 x 10° Salt water (saturated) 4.4 x 10°? 
Copper 1.68 x 10-8 Germanium 4.6-10' 
Gold 2.21 x 1078 Diamond 2.7 
Aluminum 2.65 x 10°8 Silicon 2.5 x 10° 
Tron 9.61-10°° Insulators: 

Mercury 9.58 x 1077 Water (pure) 2.5 x 10° 
Nichrome 1.00 x 10°° Wood 10°10"! 
Manganese 1.44 x 10°° Glass 10'°-10'4 
Graphite 1.4 x 107> Quartz (fused) ~ 10'° 
Solution 


(a) Equation 4.99 is +7 = e/o. But € = €o€, and from Eq. 4.52, we have e, & n 


> 


and for glass the index of refraction is typically around 1.5, so e ~ (1,5) 
8.85 x 10° '*C?/Nm’, while o = I/p ~ 10° '7Qm (Table 4.3). Then t = (2 x 
10~''y/10~!? = 20s. (But the resistivity of glass varies enormously from one 
type to another, so this answer could be off by a factor of 100 in either direction.) 


(b) For silver, p = 1.59 x 10-8 (Table 4.3), and € & &, so we = 2a x 10!° x 


8.85 x 10~'? = 0.56. Because o = I/p = 6.25 x 10’Qm > oe, the skin depth 


‘ tar, fo 2 _ 
(see Eg. 4.108) is d=, Va a ——— = 6.4x 


10°-7’m=6.4 x 10-*mm. We place plate silver at a depth of about 
0.001 mm; there is no point in making it any thicker, because the fields do not 
penetrate much beyond this anyway. 








(c) For copper, Table 4.3 gives o = 1/(1.68 x 10-°) +6 x 10”, wep = (2a x 10°) x 


(8.85 x 10!) = 885 x 10°. Because o > we, Eq. 4.106. > kx 


2m (P= 2 | aoeciagsiet = 4% 104 m= 04 mm. From 
Eq. 4.107, the propagation speed is v = 2 = 3A = Av (4 x 10-*) x 10° = 
400 m/s. In vacuum, A = c= be = 300m; 0 =c=4 x 10°m/s. (But really, in 
a good conductor the skin depth is so small, compared to the wavelength, that the 
notions of “wavelength” and “propagation speed” lose their meaning. 




















4.9.2 Reflection at a Conducting Surface 


The boundary conditions we used to analyze reflection and refraction at an interface 
between two dielectrics do not hold in the presence of free charges and currents. 
Instead, we have the more general relations (see Eq. 4.51) with addition of the term 


for surface current, K f in (iv) [1]: 
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(i) e£t—e£t=0 (iii) E!-zl =o 
1 L : lo 1 al Za + (eles) 
(ii) By — By =0 (iv) rad aay = (Ky) x (7) 


In this case o(not to be confused with conductivity) is the free surface charge K fs 
the free surface current, and ii (not to be confused with polarization of the wave) is a 
unit vector perpendicular to the surface, pointing from medium (2) into medium (1). 
For ohmic conductors, (J, p= cE), there can be no free surface current because this 
would require an infinite electric field at the boundary [1]. 

Suppose now that the xy-plane forms the boundary between a non-conducting 
linear medium (1) and a conductor (2). A monochromatic plane wave, traveling in 
the z-direction and polarized in the x-direction, approaches from the left, as shown 
earlier in Fig. 4.20: 


= ee ts a s | a 
E,(z,t) = Ee) % and B7(z,t) = — Epc ¥ (4.124) 
vi 
This incident wave gives rise to the following reflected wave [1]: 
i ~ aa Ay = 1 = i(—k,z-at) > 
Ep(z,t) = Eo,e"-X and Br(z,t)= — 5, Fone way (4.125) 
1 


Propagating back to the left in medium (1) produces a transmitted wave: 


zp LE i(k a a ko ~ fC 
Er(z,t) = fig, eiliae-o1) and Br (é2) = 2 By ellh art) 


yz (4.126) 


which is attenuated as it penetrates into the conductor. 

At z = 0 the combined wave in medium (1) must join the wave in medium (2), 
pursuant to the boundary conditions (see Eq. 4.123). Because E~ = 0 on both sides, 
boundary condition (i) yields of = 0. Because B* = 0 (ii) is automatically satisfied. 
Meanwhile, (iii) gives: 





Eo, + Eog = Eo, (4.127) 
and (iv) (with Ky = 0) says, 
: (Eo, — Eve) ke Ep, = 0 (4.128) 
0, Or — ‘: 
Myvi “yo 
or 
Eo, — Eq, = BEo, (4.129) 


where 
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pa (4.130) 
ao) 
It follows that 
7 LaF ys ~ 27\. 
Ba =(TF)en ant Bo, = (2) Eo winor 


These results are formally identical to the ones that apply at the boundary between 
non-conductors (see Eq. 4.66), but the resemblance is deceptive because B is now a 
complex number. 

For a perfect conductor (6 = 00), kz = oo (see Eq. 4.105), so B is infinite, and 


Er => —Eo, and Eo, =0 (4.132) 


In this case the wave is totally reflected, with a 180° phase shift. (That is why 
excellent conductors make good mirrors.) In practice, you paint a thin coating of 
silver onto the back of a pane of glass. The glass has nothing to do with the 
reflection; it is just there to support the silver and to keep it from tarnishing. Because 
the skin depth in silver at optical frequencies is on the order of 100 A°, you do not 
need a very thick layer. 


Example 4.10 Calculate the reflection coefficient for light at an air-to-silver inter- 
face (4; = 2 = fo, €1 = €, 6 = 6 x 10’(Q - m)"') at optical frequencies 
(o =4 x 10574). 

Bop? 
Ey, 





Solution According to Eq. 4.123,R 








iD ss re 
1-p 1-p)\ (1-p* axe: 
is (+4) (1552) where " 


a complex conjugate of # and fp = alka. ky = ala (kp + ik) (see Eqs. 4.110 and 
H2@ H2@ 











4.130). Because silver is a good conductor (o A Pi Eq. 4.111 reduces to 


— Yk W,/ Se a 4 (=e, sO p= aaa tk = ac lag + ik2). 
H2@ 
[Sto 3x10") (42x 10- 
Now let Y= MyV14/ a = Hoc \/ aa cy/ Fe = (3 * 10°) ( ea f = 


29 and R= (; Y ) (; =i) == aA = 0.93. Evidently 93% of the light is 


l+y+iy l+y—iy 














reflected. 


4.9.3 The Frequency Dependence of Permittivity 


In the preceding sections we have seen that the propagation of EM waves through 
matter is governed by three properties of the material, which we took to be constants; 
they are as follows: (1) permittivity ¢, (2) permeability, and (3) conductivity o. 
Actually, each of these parameters depends to some extent on the frequency of the 
waves you are considering. 
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Fig. 4.27 Graph for typical 1.480 | 
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Fig. 4.28 Phase and group velocity [1] 


Indeed, if the permittivity were truly constant, then the index of refraction in a 
transparent medium, n = ,/é,, also would be constant. Nonetheless, it is well known 
from optics that n is a function of wavelength. Figure 4.27 shows a graph for typical 
glass. A prism or a raindrop bends blue light more sharply than red and spreads white 
light out into a colorful rainbow. This phenomenon is called dispersion. By exten- 
sion, whenever the speed of a wave depends on its frequency, the supporting 
medium is called dispersive. Conductors, incidentally, are dispersive—see 
Eqs. 4.110 and 4.111. 

Because waves of dissimilar frequencies travel at various speeds in a dispersive 
medium, a wave form that incorporates a range of frequencies will change shape as it 
propagates. A sharply peaked wave typically flattens out, whereas each sinusoidal 
component travels at the ordinary wave (or phase) velocity, 


) 
=— 4.133 

v=" (4.133) 

The packet as a whole (the “envelope’”) travels at the so-called group velocity [9]. 


For more information, refer to Appendix G of this book. 


_ do 


== (4.134) 


Ug 


Figure 4.28 is a typical depiction of these two wave velocities. 
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Fig. 4.29 Electron 
movement [1] 


Electron 


. K spring 











Our purpose in this section is to account for the frequency dependence of ¢ in 
non-conductors, using a simplified model for the behavior of electrons in dielectrics. 
The electrons in non-conductors are bound to specific molecules [1]. The actual 
binding forces can be quite complicated, but we will picture each electron as if 
attached to the end of an imaginary spring, with force constant Kgpring (Fig. 4.29): 

Foinding —_ —KspringX — —manx (4.135) 
where x is displacement from equilibrium, m is the electron’s mass, and @po is the 
natural oscillation frequency, \/Kspring/m. Using the Taylor series and expanding the 
potential energy around the equilibrium point for sufficiently small displacement, we 
have [1]: U(x) = U(O) + xU'(0) + $x°U"(0) sb saie 

The first term is a constant, with no dynamic significance (you can always adjust 
the zero of potential energy so that U(0) = 0). The second term automatically 
vanishes because dU/dx = — F, and by the nature of an equilibrium, the force at 
that point is zero. The third term is precisely the potential energy of a spring with 
force constant Kepring = @Uldx" ly (the second derivative is positive for a point of 
stable equilibrium). If the displacements are small, the higher terms in the series can 
be neglected [1]. Meanwhile, there will presumably be some damping force on the 
electron: 


dx 
Fgamoing = is (4.136) 
The damping must be opposite in direction to the velocity and making it propor- 
tional to the velocity is the easiest way to accomplish this. The cause of damping 
does not concern us here—among other things, an oscillating charge radiates and the 
radiation siphons off energy. 
In the presence of an EM wave of frequency w, polarized in x-direction (see 
Fig. 4.29), the electron is subject to a driving force 


Fariving = QE = qE, cos (at) (4.137) 


where q is the charge of the electron and Eo is the amplitude of the wave at the point 
z, which is where the electron is situated. If we concentrate only on one point and one 
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point only so that the maximum EF occurs there at t= 0, by Newton’s second law, we 
have: 


wx 
arr = Frotal = Poinding or F damping + Fariving 
(4.138) 
Oy i ein rane 
m—, my— MOW\X = COS (@ 
de Vit 0 qEo 


Our model, then, describes the electron as a damped harmonic oscillator, driven at 
frequency @ [1]. (The assumption is that the many more massive nuclei remain 
at rest.) 
Equation 4.138 is easier to handle if we regard it as the real part of a complex 
equation: 
aR WE pe 4 yet 4.139 
oo oe 0e (4.139) 


In a steady state the system oscillates at the driving frequency: 
x(t) =Xpe7™ (4.140) 
Inserting this into Eq. 4.139, we obtain 


a q/m 
xo = 5} 
Qo 


0 yo 


The dipole moment is the real part of the following: 
Eye” (4.142) 


The imaginary term in the denominator means that p is out of phase with E—lagging 
behind by angle tan ~! [yo / (w5 a w)| , which is very small when @ < @o and rises 
to z when @ > @p. 

In general, inversely situated electrons within a given molecule experience 
various natural frequencies and damping coefficients. Let us say there are f; electrons 
with frequency w; and damping y; in each molecule. If there are N molecules per unit 
volume, the polarization P is given by the real part of 


~  Nq f; = 
P= aft E 4.143 
m (>: wy —@ — =a ( ) 


This applies directly to the case of a dilute gas; for denser materials the theory is 
modified slightly, in accordance with the Clausius-Mossotti equation. Note that we 
should not confuse the “polarization” of a medium, P, with the “polarization” of a 
wave—same word, but two completely unrelated meanings. 
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N ow we define the electric susceptibility as the proportionality constant between 
Pand E (specifically, P= EN E) [1]. In the present case P is not proportional to E 
(this is not, strictly speaking, a linear medium) because of the difference in phase. 
The complex polarization P, however, is proportional to the complex field E, and this 
suggests that we introduced a complex susceptibility, y ,: 


P = eny,E (4.144) 
From all these we conclude that physical polarization is the real part of P, just as 
the physical field is the real part of E. Specifically, the proportionality between D 


and E is the complex permittivity, € = eo(1 + var and the complex dielectric 
constant in this model is: 





: N@’ ie 
ae j 4.145 
: me a rere ( ) 


Ordinarily, the imaginary term is negligible; however, when @ is very close to one 
of the resonant frequencies (q@;), it plays an important role, as we will see. In a 
dispersive medium the wave equation for a given frequency reads: 


WE = Buy oe — (4.146) 


This provides a plane wave solution, as before: 


F(z, 1) = Eel) (4.147) 
with the complex wave number 
k = Jeo (4.148) 
Writing k in terms of its real and imaginary parts, 
k=k+ix (4.149) 
Equation 4.147 becomes 
E(z,t) = Eye “ello (4.150) 


Evidently the wave is attenuated, which is not surprising, because the damping 
absorbs energy [1]. Because the intensity is proportional to E” and thus to ye **“, the 
quantity 


a=2k (4.151) 
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is called the absorption coefficient. Meanwhile, the wave velocity is w/k, and the 
index of refraction is: 


n=— (4.152) 
@ 


Nevertheless, in the present case k and x have nothing to do with conductivity; 
rather, they are determined by the parameters of the damped harmonic oscillator. For 
gases the second term in Eq. 4.145 is small, and we can approximate the square root 
(see Eq. 4.148) by the first term in the binomial expansion, V1 + ¢ = 1 + 5é. Then, 














> oO —_ oO Ng’ i 
k= r& 1 4.153 
ee | ne 2 =a ee ( ) 
J J J 
or 
2 8 
n=—e14 yo (4.154) 
ad ee a (03 - «”) — iva 
and 
Nee Yi 
a=2e214-2"° 5° fy (4.155) 


Griffiths plotted the index of refraction and the absorption coefficient in the 
vicinity of one of the resonances [1]. Most of the time the index of refraction rises 
gradually with increasing frequency, consistent with our experience from optics (see 
Fig. 4.26). In the immediate neighborhood of a resonance, however, the index of 
refraction drops sharply. Because this behavior is atypical, it is called anomalous 
dispersion. Notice that the region of anomalous dispersion (@, < @ < @» in the 
figure) coincides with the region of maximum absorption; in fact, the material may 
be practically opaque in this frequency range. The reason is that we are now driving 
the electrons at their “favorite” frequency; the amplitude of their oscillation is 
relatively large, and a correspondingly large amount of energy is dissipated by the 
damping mechanism. 

In Fig. 4.30 n runs below | above the resonance, suggesting that the wave speed 
exceeds c. This is no cause for alarm because energy does not travel at the wave 
velocity but rather at the group velocity (see Example 4.11). Moreover, the graph 
does not include the contributions of other terms in the sum, which add a relatively 
constant “background” that, in some cases, keeps n > 1 on both sides of the 
resonance. 

If you agree to stay away from the resonances, the damping can be ignored, and 
the formula for the index of refraction simplifies: 
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Fig. 4.30 Plot of the index 
of refraction and absorption 
coefficient near of one of the 
resonances [1] 








Ng fi 

n= 1+— =F 4.156 

Meo » o; — @ ( ) 

For most substances the natural frequencies, w,, are scattered all over the spec- 

trum in a rather chaotic fashion. But for transparent materials, the nearest significant 
resonances typically lie in the ultraviolet, so w < @;. In that case, 


1 1 wo 1 a 
w—o we o = er 
ij ij j j 





nN 


and Eq. 4.148 takes the form 





Ng fj 2 Ne fi 
pe (ee ee ae ee ee 4.157 
‘“ meo » oo; 7 2meEo » ws ( ) 
Or, in terms of the wavelength in vacuum, (J = 2zc/o): 
B 
n=1+A Bag (4.158) 


This is known as Cauchy’s formula; the constant A is called the coefficient of 
refraction and B is called the coefficient of dispersion. Cauchy’s equation applies 
reasonably well to most gases in the optical region. 

What Griffiths described is certainly not the complete story of dispersion in 
non-conducting media [1]. Nevertheless, it does indicate how the damped harmonic 
motion of electrons can account for the frequency dependence of the index of 
refraction, and it explains why n is ordinarily a slowly increasing function of @, 
with occasional “anomalous” regions where it precipitously drops. 
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Example 4.11 Assuming negligible damping (y; = 0), calculate the group velocity 
(v¢ = da) of the waves described by Eqs. 4.150 and 4.153. Show that v, < c, even 
when v > c. 
Solution 
Nq d 1 

cd fi and v, = Be 


dk (dk/da) 








dk Cc 2meo = w) (03 — 
2 2 
1 Ne xr ( ;+@ ) 
= 2meo j 2 2 2 
(05 - «?) 
(0? + w) 





Because the second term in square brackets is positive, it follows that v, < c, 
f 
2 
2 


(i 





a) Nq@ ‘ 
whereas v = ——< 1+ > is greater than c or less than c, 


2meo 


depending on . 


4.10 Electromagnetic Waves in Conductors 


In previous discussions of EM waves in matter, we have been assuming there are no 
free charges (Pfree = py) Or free currents (J free = J, Ar the only charges present were 
the bound charges in the dielectric. When dealing with conductors, there are plenty 
of electrons that are free to move about when an external electric field is applied. 
There are still bound charges in a conductor (atoms in a conductor generally only 
contribute one or two conduction electrons each), but we will find that (in a good 
conductor) the interaction of the fields with the free charges and currents will 
dominate everything else. 
In an ohmic materials that obeys Ohm’s Law, we can write: 


IG) =o. Eg) (4.159) 


where o, is the conductivity of the conducting materials (not to be mistaken for the 
surface charge density!) and o, = 1/p,, which we will take to be uniform inside the 
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conductor (i.e., 9. = resistivity of the metal conductor with a dimension of Ohm- 
meter). For any free charge that might build up in the conductor because of these 
currents, we know from charge conservation: 


9-H) =- 22 (4.160) 


Thus, inside such a conductor, we can assume that the linear/homogeneous/ 
isotropic conducting medium has electric permittivity e and magnetic permeability 
uu. Combine this with Gauss’s Law for electric fields in matter, 


vo (cE (7, t)) =—py, to find how quickly the free charge distribution would 
dissipate. Bear in mind that Maxwell’s equations always are true, in matter or not; 
splitting things up into free and bound charges (or for that matter current) can be 


convenient for solving for the fields, but we do not have to do things this way. Now 
substituting these relations into Eq. 4.152, we can write: 


Pr _ %. (eB(7,1)) = -% [(e8(7,))] =o, (4.161a) 
or 
ws <p; (4.161b) 


The solution to this first-order differential equation is: 


p(t) = p(0)exp(—o.t/e) = p(7,t = 0)e1/™™ (4.162) 


This is a damped exponential type of function. This also tells us that any free 
charge that might build up in a conductor will dissipate (because of the mutual 
repulsion of the charges) with time constant 7, = é/o,. The “time constant” repre- 
sents the amount of time for the initial density to reduce to 1/e(~37%) of its original 
value. For a perfect conductor (6, — 00), the time constant goes to zero, meaning the 
charge density instantly dissipates. For a “good” conductor (e.g., copper), where 
the conductivity is typically on the order of 10° (Q m)~’, the time constant is 
somewhere around 10~'° s. This is actually much smaller than the typical time 
between collisions of the electrons with the atoms making up the conductor, which is 
t.~107 14. 

For frequencies higher than 1/r,, Ohm’s Law starts to break down, so 7, is the 
important time constant here. Therefore, the assumption of dealing with the linear/ 
homogeneous/isotropic conducting medium is a valid assumption for the rest of this 
discussion. Then, let us assume we are working with good conductors with frequen- 
cies that are at or below the optical range (10'° Hz), which is pushing the use of 
Ohm’s Law, but not so much that the results are no good. In other words, we should 
have @ < o-./e or ew@/o, < 1. It turns out that this assumption will make the 
calculations much easier as we go along. 
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Note that going for forward with further analysis of EM wave through conduction 
we write the symbol o, for conductivity of conducting materials as o,; now that we 
know that is the case it is not indeed the surface charge density. Also, we start using 
vector notation rather than bold characters. 

Now again assuming a conductor of linear, homogeneous, and isotropic 
conducting material that has electric permittivity e and magnetic permeability yp, 
we substitute py = 0, which means free charges dissipate rather quickly and J (7, t) 
=oE ( fr t) This also means that J ( r, t) # 0, then we can write Maxwell’s equations 
inside such a conductor are as follows: 





(1) o-B(,) 2 

(2) V- B(7,1) =0 (4.163a) 
a3 OB(?, 

(3) Vx E(#,1) =- e ‘ 


and using the Ohm’s Law relation of a (7, t) =oE (F, i); we have the fourth of 
Maxwell equation as: 
OE(?,t : OE(F,t 
2B) _ 50H) 4 pe Ee 
Ot 
Because electric charge is always conserved, the continuity equation inside the 
conductor ends up with a homogeneous differential equation of first order as 


Eq. 4.161b, where the solution is provided in Eq. 4.162 and depicted in Fig. 4.31, 
which is a damped exponential curve. 


(4) V x B(#,t) =p (7,1) + pe (4.163b) 


Fig. 4.31 Sketch of free Phirecl?, t) 
density dissipated vs. 4 
characteristic time (solution 

to Eq. 4.162) 


l/e}-- 
of initial value ' 
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Example 4.12 Calculate the charge relaxation time for copper. 


Solution Assume the following data for copper: 





Pou = 1/ocu = 1.68 x 10°Q—m 6cu = 1/Pey = 5.95 x 10’Siemens/m 





If we assume €cy © 3€) = 3 x 8.85 x 10°-°F/m for copper metal, then: 


TE (€cu/ocu) = PouEcu = 4.5 x 10~°s 


However, knowing that the characteristic/mean collision time of free electros in 
pure copper is 76°! ~ 260" /ve" 4, where 20” ~ 3.9 x 10-8m = mean free path 
(between successive collisions) in pure copper, and Sea ~ V/3kgT/m,. ~ 12 x 
10°m/s; thus, we obtain 72" ~ 3.2 x 1071°s. 

Consequently , we see that the calculated charge relaxation time in pure copper, 
a ~ 4.5 x 107!°s, is much less than the calculated collision in pure copper, 
ty ~ 3.2 x 10's. Furthermore, the experimentally measured charge relaxation 
time in pure copper is TS (experimenta) ~ 4.0 x 10°!°s, which is roughly a 
fivefold magnitude greater than the calculated charge relaxation time, 
to AS 10s: 


The problem here is that the macroscopic Ohm’s Law is simply out of its range of 
validity on such short time scales! Two additional facts are that both, and C 4 6, are 
frequency-dependent quantities (i.e, ¢€ = e(@) and o = o(@)), which becomes 
increasingly important at the higher frequencies (f = 22/@~1/tye1ax) associated 
with short time-scale, transient-type phenomena! 

So, in reality, if we are willing to wait even a short time (e.g., At~1 ps = 10~'7s), 
any initial free charge density, Pcl? t= 0), accumulated inside the conductor at 
t = 0 will have dissipated away (i.e., damped out), and from that time onwards, Prog 
(F, t) = 0. can be safely assumed. 

Thus, after many charge-relaxation time constants (e.g., 20T;e1a, << At ~ 1 
ps = 10~'?s), Maxwell’s equations for a conductor become (with Pree (7, t > At) 
= 0 from then onward); that is, new sets of Maxwell’s equations for a charged- 
equilibrated conductors: 


(1) V- E(7,2) =0 
(2) V- B(?,1) =0 
(3) ox E(7,) = - AD 
ee 7 OE (7. we r 
(4) V x B(F,t) = pd (7,t) +e Ao = ob 7) 4A) 
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Now because these equations are different from the previous derivation(s) of 
monochromatic plane EM waves, propagating in free space/vacuum and/or in linear/ 
homogeneous/isotropic non-conducting materials (only Eq. (4) has changed), we 
again derive the wave equations for Eand B from scratch. As before, we apply acurl 
operation to Eqs. (3) and (4) and using the following vector identity, we get: 


> a 


Vx(Vx E)=(V-E)V—-(V-V)E+(E-V)V—(V-V)E (4.165) 



































Ox (Fx A) =-L£ (9 xB) 
= V( y BE)" —-WE= 5 nod ve SE) 
S545 => Or OE 
_—W2FHE 
Nene On | Mo Or 
(4.166) 
V x (V x B) = plo(¥ x B)| $en(9 x E) 
cee do SOR... 2058 
= VV B) — V°B = po> — peaz 
2 OB OB 
= V'B = pes + po 
Again, we can write: 
2A(2> > 
V7E(7,t) aye ae + po AB!) and 
a He (4.167) 
Us O° B(7,t OB(7,t 
2 > _ ? aL 5) 
V B(?,t) He 5 HM Fr 


Note that these four-dimensional wave equations for E and B in a conductor have an 
additional term that has a single time derivative, which is analogous to a velocity- 
dependent damping term (e.g., for a mechanical harmonic oscillator). 

The general solution(s) to the preceding wave equations are usually in the form of 
an oscillatory function, *, a damping term (i.e., a decaying exponential) in the 
direction of the propagation of the EM wave; for example, complex plane-wave 
type solutions for E and B associated with the preceding wave equation(s) are of the 
general form: 


E(z,t) = Ege") (4.168a) 


and 
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aa ee Ve k\~ x 
B(z,1) = Boel) = (=)i EGH= 


@ 


I~ 3s 
—k x E(z,t) (4.168b) 
wo 

With a frequency-dependent complex wave number, k(w) = k(@) + ix(w), where 
k(w) = Re(k(w)) and x(w) = Im(k(@)); thus, the corresponding complex vector 
in the positive +2 -direction is k(w) =k(w)k =k(w)Z —that is, k(w) = 
[k(@) + ik(@)|z. 

We plugged E(z,t) = Ege!(&-o1) and B(z,t) = Bye! (eo) into their respective 
wave equations before, and obtained from each the same/identical characteristic 
equation (i.e., dispersion relation) between complex k(@) and @; thus, we get the 
following relation: 





k?(w) = pew? + inow (4.169) 


Thus, because k(w) = k(@) + ix(@), then: 








k?(@) = [k(@) + ix(w)|° = °(@) — K?(w) + 2ik(w)K(@) = pew + iuow 
(4.170) 








If we temporarily suppress the w-dependence of complex k(a), this relation 
becomes as: 
k?(w) = (k + ix)? =? — x + 2ike = pew? + ipow (4.171) 


We can solve this relation to determine (w) = Re(k(w)) and k(@) = Im(k(@)) as 
follows. 

First, separate this relation into two relations by separating its real and imaginary 
parts: 


k—K = pew and 2ikk = inow 2kk = pow (4.172) 
Now we have two separate independent equations, k* — x? = pew” and 2kx = pow, 


and we have two unknowns (i.e., k and x). Therefore , solving these two equations 
simultaneously, we find the following results: 


1 
K= zHow/k 


2 2 _ 72 1 : 2_ If! : 2 
ke —K =k zhow/k =k gag OO) ee 


Then multiply by k and rearrange the terms to obtain the following relation: 


(4.173) 


1 2 
kt — (ew?) kK? — (5 pow =0 (4.174) 
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To solve this equation, let x = Kk, a=1,b= (uew"), andc = — (Su0w)”, then 
Eq. 4.170 reduces to the form of a quadratic equation of ax’ + bx +c = 0 with 
solution roots of 


E 1 ; ae ae a 
a a a 4( — 
5 te) 5 + (eo) | oe + (Jucw) 


1 (u°0"w") 
k* = ~(uom") |1 = 4/1+4 
(uow”) 4(u2e2o") 


7 (uo) r 24s (2) ] 


Now we can see that on physical grounds (k° > 0), we must select the + sign, thus: 


14 +) 



































Pas 


=} (uo«?) 





and thus: 


k= VP =0 at bal 14 2) ne c + 4/14 | - (4.175) 


Having solved for k (or equivalently k°), we can use either of the original two 
relations to solve for x (e. g. k? — x* = pew”), then: 











K =k — pew = = (uew’) 








1 
2 


Thus, we obtain: 


(4.177) 





Note that the imaginary part of k, x(@) = Im(k(w)), results in an exponential 
attenuation and damping of the monochromatic plane EM wave, increasing z as: 


ei Eye-“el- for Electric Field 


ty 


(4.178) 


joel) 


a . 11> 3 = 
(z, t) = Boeel(-@) — kx Egee'(*-©) for Magnetic Filed 
@ 
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These solutions (Eq. 4.178) satisfy the preceding wave equations for any choice 
of E 0: 

The characteristic distance over which E and B are attenuated, and reduced to 
1/e = e~' = 0.3679 of their initial value at z = 0, is known as the skin depth and is 
shown as 6¢kin depth = Osc = 1/k(@) and the SI unit has a dimension of the meter: 





(z = bg, t) = Ee !eil—-o0) 
(Z = bcc, t) = Boe7 eile) 





jas 
Don Pye 











(4.179) 


The real part of k (i.e. k(w) = Re(k(w)) ) determines spatial wavelength A(@), 
propagation speed v(@) of the monochromatic EM plane wave (see Fig. 4.25) in the 
conductor, and also the index of refraction: 














ia) = 2m 2a 
1 Ka) ‘Re(k(o)) 
v(@) = Ko) = He(k()) (4.180) 
ie = @ —ck(@) | cMRe(k(w)) 
v(@) 7) @ 


7 The preceding plane-wave solutions satisfy the wave equations for any choice of 
Eo. As we also have seen, it can similarly be shown here that Maxwell’s equations 
(1) and (2), iV E=Oand V-B= 0}, rule out the presence of any (longitudinal) 
z-components for E and B for EM waves propagating in the +Z-direction: > Eand B 
are purely transverse waves. 

If we consider, for example, linearly polarized monochromatic plane EM waves 
propagating in the + -direction in a conducting medium (e.g., E(z,t) = Epe~* 
el(k-o)2), then: 


2 Ts. 3 k = ; k+ik\~ —. on 
B(z,t) = ok x E(z,t) = (5) Ege Melle o5 = (i *) Ege el2-o)5 


1 
ty 
o 
Kk 
joel 


(z,t)LZ (+2 = propagation direction) 
(4.181) 


The complex wave number is k =k+ ik = Ke'*, where K= |k| = =VP+ 
and @, = tan “(k/k) (Fig. 4.32). Then, we see that: E(z, t)= Eye- Kei(k—ot)F has Eo 
= Epe’** and that: 
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Fig. 4.32 In the complex k 
-plane 





Blz t) = Boe “ell OF = Xe ell-oy 


~ Kei* is 
Ey = Eoe’£. Thus, we see that: 
oy) 


af ge 
OT rpelletoe) 
@ 





LPy 
Ke" Eye® = K Epel(Get os) = 
o 
That is, inside a conductor E and B are no longer in phase with each other. Phases of 


E and Bare dn = de + b | We also see that 6g _ ¢ = 6g — 6g = &; magnetic field 


lags behind the electric field. We also see that: 


B 2 1 
2a En |e ny re 
Eo @ F200) Cc 


The real physical E and B fields associated with linearly polarized monochromatic 
plane EM waves propagating in a conducting medium are exponentially damped: 














Ca Re( Ete, ‘)) = Ege cos (kz — wt + 5g)x (4.182a) 
B(z, t) = Re(B(z, 1) = Boe cos (kz — wt + 6g)5 (4.182b) 
= Bo cos (kz — wt + {6 + o})¥ . 
Bo K(a) 2] 
Pa aaa an 
= Pree (=) | (4.183a) 





1/2 
K(@) = e(o)| = /K (0) + (@) = of 1+ (=) ] (4.183b) 
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The definition of the skin depth in a conductor is: 





Distance over which, 
= 7a =| the Eand Bfields fall to 
€| o 2 by ek Baie a 
or/F | Mog (S) | their initial values 








5(@) cin depth = k(@) 




















(4.185) 


Example 4.13 What is the skin depth of fine silver at a microwave frequency of 
10!° Hz, which is the common microwave region; assume the silver has a 
conductivity of g = 10'° Hz. 


Solution The skin depth is presented by the following equation: 


= 2 2 7 , 
oS le 7 I x 10'°) (4x x 10) (3 y 10°) =92x 10cm 


Thus, at microwave frequencies the skin depth in silver is very small, and, 
consequently, the difference in performance between a pure silver component and 
a silver-plated brass component would be expected to be negligible. 








Example 4.14 For a seawater case we calculate the frequency at which the skin 
depth is 1 m. For seawater use “ = fy and g © 4.3 S/m. 


Solution The expression for the frequency corresponding to given skin depth 6 is: 


2) 2 3.70 x 10° 


—1 
— = — S 
OF eye 4.3.x 4a x 1078 ge 





which yields f= 58 x 10° HZ, or frequency of 60 kHz for a skin depth of 1 m. Ifa 
submarine is equipped with a very sensitive receiver, and if a very powerful 
transmitter is used, it is possible to communicate with a submerged submarine. 
However, a low radio frequency must be used, and even then, an extremely 
severe attenuation of the signal occurs. At five skin depths (5 m in the case calculated 
here), only 1% of the initial electric field remains and only 0.01% of the incident 
power. 


412 4 The Fundamental of Electrodynamics 
References 
1. D. Griffiths, Introduction to Electrodynamics, 3rd edn. (Prentice Hall, New York, 1999) 


2. R.M. Lambert et al., J. Phys. E: Sci. Instrum., 2799-2801. https://doi.org/10.1088/0022-3735/2/ 


ONAUAA 


9/311 

.J.R. Reitz, FJ. Milford, R.W. Christy, Foundations of Electromagnetic Theory, 3rd edn. 
(Addison Wesley, Reading, 1979). Section 17-5 

. L. Eyges, The Classical Electromagnetic Field (Dover Publication, New York, 1972) 

M.B. James, D.J. Griffiths, Am. J. Physics 60, 309 (1992) 

.N. Ashby, Am. J. Phys. 43, 553 (1975) 

H.C. Ohanian, Am. J. Phys. 51, 1020 (1983) 

. W. M. Hayes (ed.), Handbook of Chemistry and Physics, 78th edn. (CRC Press, Inc., Boca 
Raton, 1997) 

.F.A. Jenkins, Fundamentals of Optics, 4th edn. (McGraw-Hill Science/Engineering/Math, 
New York, 2001) 


